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(SADOI or DOI) 

Closed 

fixed, 42, kj>, kk relative 
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22 + 2n milliseconds 

Adjacent to temporary storage 3^ 

Adjacent to routine 17 

Within routine 8 

The function, f (x), must be given in a sequence of n values 

(at least 15 ) for equal increments, £x s of the argument, x, 

and must be such that -1 < f < 1 for all i> 1 < i < 15 < n 

where f. m f(x. ) Furthermore the range x, to x must be 
i l J. n 



such that x - x, 
1 n 1 ' 



|f 1 l<2. 



f f (x) cjx is approximated by the formula* 

1/2 j f M <** = ** f a i ( ? r i tf n ) + a 2 ( f 2 +f n-l J + 
x l 

a 3 (f 3 +f n-2> + WW + WW + 

"WW + WW + & ( t 6*'" +t a . 7 )i 
where the a's are such that f(x) is approximated by a 6th 
degree polynomial. 

Entry to the subroutine is made vi% orders 
J-1 L5 n-14 L5 A 
4 L5 B L5 J 
j+1 26 m 
where n> 15 is the number of tabulated function values, 
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A is the address of the first of these, f , B is the address 



of the increment, Ax = x ^ - x., and m is the address of 
the first word of the integration routine. Control is rei 
to the right hand side of j+1 with l/2 f n f fojdx in R . 



*G. Bland and I. Rhodes, "Seven-point Lagrangian Integration 
Formulas", Math. Table MS 25, U.S. Department of Commerce. 19^9 
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